Haldane Gap and Fractional Oscillations in Gated Josephson Arrays 
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An analogy between the twisted quantum xxz model and a gated Josephson junction array is used 
to predict sharp structure in the critical currents versus gate voltage, and fractional AC Josephson 
frequencies. We prove selection rules for level crossings which imply fractional periodicities of ground 
states with varying Aharonov-Bohm flux. Extrapolated numerical diagonalization on ladders, find 
a Haldane gap at moderate easy-plane anisotropy, with vanishing superfluid stiffness. Physical 
parameters for experimental realization of these novel effects are proposed. 
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Quantum phase fluctuations in superconductors can 
drive zero temperature superconductor to insulator tran- 
sitions fl, as observed for example in disordered films 
P). Their effects can be enhanced and studied in de- 
tail using a weakly coupled, low capacitance Josephson 
junction array (JJA) ffl]. 

Theory of quantum phase fluctuations has used models 
of interacting bosons on a lattice |5| , and quantum dy- 
namics of vortices ||. The latter approaches have pro- 
posed collective phases such as vortex Bose condensa- 
tion Q (for the insulator), and fractional Quantum Hall 
phases (for JJAs in magnetic field ||). 

Lattice bosons map onto effective models of quantum 
spins. A popular approximation to the phase diagram is 
mean field theory on the classical (large S) spin model 
In the strongly quantum regime, the same mapping 
relates the Mott insulator (integer bosons per site) and 
the quantum disordered antiferromagnet p]. 

In this paper we explore the quantum magnetism anal- 
ogy further. We focus our attention to the effects of a 
periodic lattice on superconductivity. We study the quan- 
tum xxz model with twisted boundary conditions (i.e. an 
Aharonov Bohm (AB) flux) both numerically and analyt- 
ically. The many body spectrum, vortex tunneling rates, 
and superfluid stiffness are computed for different lattice 
dimensions and magnetization (Cooper pair density). We 
prove general selection rules for symmetry protected level 
crossings. This rule imposes fractional periodicities of 
the ground state as a function of AB flux, and is closely 
related to the "fractionally quantized phases" found by 
Oshikawa, Yamanaka and Affleck (OYA) in magnetized 
Heisenberg chains JlO[ . 

The following effects may be observed in JJAs of di- 
mensions L x x Ly. 

(i) Haldane gap. Our numerical results for the xxz model 
on a two leg ladder find a Haldane gap in a sizeable 
regime of the easy plane anisotropy parameter. The low- 
est gap remains finite, while the superfluid stiffness de- 
cays exponentially for L x — > oo. The Haldane phase is 
characterized by a suppressed critical current (relative to 
its classical value), and a high AC Josephson frequency 



Jq = 2eV/h, where the classical Josephson frequency of 
the array is f c i = 2eV/{hL x ). 

(ii) Fractional Oscillations. At fractional Cooper pair 
densities n = p/(qL y ), p,q integers, selection rules de- 
rived below produce sharp dips in the critical current I cr 
versus gate potential. In these states, AC Josephson os- 
cillations appear at subharmonic frequencies /q = Jq jq. 

We conclude by proposing physical parameters for ex- 
periments. 

The short range Bose Hubbard model is given by 
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where b\ creates a boson (Cooper pair p|) at site i on 
a square lattice with nearest neighbor bonds {ij}, and 
n-i = b\bi. The lattice is placed on a cylinder penetrated 
by an AB flux introduced via the gauge phases % = 
6j t i +x 2ir<j>/L x , where <f> = $/$o, and $ = h/(2ec). The 
supercurrent in the x direction is given by I s — i ( & } . 

At large U » J, V one can keep the two lowest energy 
Fock states at every site, say |n$) and |flj + l), and project 
out all other occupations. In the projected subspace, 
b\,bi,ni — Hi are replaced by spin half operators S^~, S~_, 
and Sf + 5 respectively. This transformation maps 
onto the quantum S = \ xxz Model Il2] 
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where the Ising coupling is J z = V/A, and we limit our 
discussion to easy plane anisotropy J z < J, in order to 
avoid the charge density wave phases Q, A pure gauge 
transformation <f> — > <fi + 1 on ([!]) or (||) leaves their spec- 
trum invariant. 

It is instructive to consider the classical (large S) 
ground state energies of (||) which are adiabatically con- 
nected to the ground states at <f> = i, i = 0, 1, . . ., 
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E l cl (cj)) = JL y L x 5n(l 
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where Sn = n — h. At <j) — q/2, q integer, pairs of clas- 
sical ground states of oppositely directed supercurrents 
become degenerate, as depicted in Fig. dll). 
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FIG. 2. Haldane gaps. Magnon excitation energies for 
two leg Josephson ladders for different anisotropy J z / J. Here 
J = 1 and i? m is the lowest eigenenergy of magnetization m 
at flux $ = $o/2. £ is the stiffness correlation length given 
by finite size scaling of ladder lengths upto L x = 12. 



FIG. I. Schematic adiabatic ground state energies as a 
function of Aharonov Bohm flux E($), where Eq = E(0). 
Thin lines: classical energies E % cl (|^). Thick lines: quantum 
adiabats of periods <?4>o, for 5 = 1,2, 3. Notice the level cross- 
ings for q 7^ I , which are protected by the selection rules (Q). 

Tunneling paths between the two classical ground 
states can be constructed as histories of vortices travers- 
ing the lattice in the y direction, or nucleation and sep- 
aration of a vortex-antivortex pairs upto the two edges. 
The tunneling matrix elements, unless prohibited by se- 
lection rules, open minigaps at the avoided level cross- 
ings. This allows their precise computation from the 
(many body) eigenenergies as a function of AB flux 

The Hamiltonian (0) was diagonalized using the Lanc- 
zos algorithm using lattice momentum and total magne- 
tization to block-diagonalize the matrix. For the two leg 
ladder at S% ot = 0, we find a regime of J z < J where 
the minigaps at <p = 1/2 remain finite as L x increases, 
which indicates a Haldane gap phase in the thermody- 
namic limit. This phase has been previously established 
for isotropic integer spin chains ]l7| , and half odd integer 
ladders and chains at finite magnetic fields (l(J. 

To establish this phase we plot in Fig. |^, the magnon 
gap |L9| at <f> = 1/2 as a function of 1/L X , and extrapolate 
the results to 1/L X — ► 0. However, for the extrapoltion 
to be justified, we must be certain that we have reached 
the asymptotic L x >> £ regime where £ is the correla- 
tion length. £ was calculated from the supcrfluid stiffness 
K = d 2 E/d<fi 2 , which was found to fit to an exponen- 
tial K(L X ) oc exp(— L x /£). We find that in the regime 



J z /J x € (0.5, 1), the correlation length is £ E (4.26, 2.13) 
which is safely smaller than the larger system dimensions. 
For the pure xy model at J z — 0, however, £ reaches our 
largest system size. In this regime therefore, an extrapo- 
lated finite gap at 1/L X = is not credible. 

An easy-plane Haldane phase is explained as follows. 
The path integral of an even leg ladder of S = 1/2 spins, 
can be mapped onto a classical partition function of an 
0(2) relativistic field theory in two dimensions |l7| . |l3| . 
Its temperature T 2d scales asymptotically as ~ (SLy)" 1 
pl| . This suggests that below a certain spin size and 
ladder width, a disordered phase with exponentially de- 
caying correlations is possible, which translates into a 
finite gap for excitations and vanishing stiffness for long 
ladders. Above a critical width L y > £^ T , correlations 
should decay as a power-law with a finite (Id) superfluid 
density 4ir 2 p s = lim^-^ (L X K ) > 0. 

Selection rules for avoided level crossings. Vortex tun- 
neling is enabled by the lattice since it breaks continuos 
translational symmetry. However, the remaining discrete 
translational symmetry imposes selection rules which are 
given by the following theorem. 

Theorem: For the Hamiltonian (^), at (j) = q/2, for in- 
teger q, any eigenstate \Sf ot ,k x ,a) where k x is the lat- 
tice momentum in the x direction, and S% ot is the total 
magnetization, is at least two-fold degenerate unless the 
following condition is satisfied: 

^ + — + SL y = p/q (p, q integers) (4) 

L x 7T 

The theorem is similar to Lieb, Shultz and Mattis (LSM) 
theorem |2(J for half odd integer spin chains, and its 
extension to finite magnetizations by OYA pfj[. Here, 
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however, we prove exact degeneracies of the twisted xxz 
model on finite lattices, while previous work concerned 
gaplessness in the thermodynamic limit at zero external 
gauge field. 

Before providing the proof, let us review three impor- 
tant classes to which the theorem for S = 1/2 applies. 

1. Odd ladders, q = 1, k x = 0, S* ot = 0. The selection 
rules (^) cannot be satisfied, implying exact ground 
state degeneracy at = 1/2. This is closely related 
to the existence of gapless excitations in the ther- 
modynamic limit of short range half odd integer 
spin chains p0| . 

2. Even ladders with integer magnetization per rung. 
The selection rule is obeyed for q = 1 , which implies 
a minigap at the first avoided crossing of the ground 
states. If this minigap survives the L x — > oo limit, 
the system is in the Haldane gap phase. 

3. Even ladders with rational magnetization per rung. 
The selection rule is obeyed only for some q > 1. 
This gives rise to a fractional AB periodicity of the 
ground state. 

Proof: The twist operator is defined as 



O(0)^exp (-i^Y^S* (r)a^ . 



(5) 



In addition, an x-inversion operator I x is defined 
^S^ylx = S% y . For any state Vo = \S% ot , k x , a), we 
define the "g-conjugated" state tp q as 



\1> q ) = d(-q)I x \1>o) 
Lemma (degeneracy): 

(1> q \H(q/2)\il> g ) = (ik\H(g/2)\1>o) 



(6) 



(7) 



The Lemma is proven by a direct substitution of (g) in 
(0) noting that 0(<f>) is the explicit gauge transformation 
on H, 



H (</>)= 6 (-<!>) H(0)6(<l>) 



(8) 



and using the identities I x O((f))I x = 0(—(p), and 
[fl"(0),J a ] = 0. 

The theorem is proved by showing that tpQ will be 
transformed by "q-conjugation" into an orthogonal state 
(V'glV'o) = unless the selection rule @ is obeyed. The 
unit lattice translation in the x direction is T x . We make 
use of the two identitites 

T x 6{-q)T- x = exp (i2nqSf ot /L x + i2nqSL y ) 0(-q), 
I x T x I x =T x l . (9) 

The lattice momentum of ip q is given by 



T x \ip q 



exp 



l ~i—S Z tot + i^qSL y + ik x ) \<ip q ) (10) 



It follows that (ipq\ipo) — unless the momentum dif- 
ference Sk x = 2ir(qSt ot / L x — qSL y ) + 2k x is an integer 
multiple of 2-7T, which proves selection rules (Q). QED. 

Translating back into the boson language, we consider 
the ground state at <f> = 0, with k x = 0, and excess 
Bose density 5n = The selection rule implies that 

as <fi is increased, the adiabatic ground state passes q — 
1 exact level crossings before reaching the first minigap 
allowed by (||). Hence it is clear that the ground state 
adiabatic periodicity in AB flux is q^o- Fig. [I] depicts 
the ground state evolution for the cases q = 1, 2, 3. The 
critical current is bounded by 



Icr(q) < {J/K)L y 5n(l 



Sn)sin(-q) (11) 

-L/ x 



which holds upto q — L x /2 where the bound coincides 
with the classical critical current lf r = (J/fr)L y 5n(l — 
Sn). 
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FIG. 3. Schematic diagram of critical currents versus 
Cooper pair density n for an array of size L X =2Q and L y =2. 
Rational numbers p/q = n q L y , which label the dips at n q , 
are indicated. The classical critical current is depicted by a 
thin line. Above the critical currents at n q , fractional AC 
Josephson frequencies are expected at 2eV/(hq). Spatial dis- 
order and finite temperature are introduced qualitatively by 
broadening of the dips. 

In Fig. U the schematic structure of the critical current 
is plotted against the Cooper pairs density for a Joseph- 
son ladder with L x — 20. Notice the sharp dips in the 
critical current at rational densities which obey the se- 
lection rule. For weak easy plane anisotropy, these min- 
ima are expected to vanish in the thermodynamic limit, 
reflecting vanishing stiffness and gapped excitations at 
these points. 

The classical Josephson frequency of an array of length 
L x array is f c i — 2eV/(hL x ) ) where V is the total voltage 
drop in the x direction. However at rational densities 
n q = p/(qL y ), for bias current slightly above I cr (n q ), an 
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AC Josephson effect should be observed with frequency 
f q = 2eV/(qh). This could be pictured as current os- 
cillations caused by moving on the adiabatic curves of 
Fig. 0. Alternatively, this effect could perhaps be better 
detected as fractional Shapiro steps in an external high 
frequency electromagnetic field p2| . 

Experimental realization. One of the experimental se- 
tups described in Ref. [Q, has individual gate voltage 
probes which control the Cooper pair density at each is- 
land to high accuracy. The short range Bose Hubbard 
model for this type of JJA can be justified if the ra- 
tio of inter-island capacitance to gate capacitance obeys 
e = C/Cq << 1. This translates to onsite and intersite 
interactions of Eq. (Q) given by U = 2e 2 (l — 4e)/Co, and 
V = e4e 2 /Co respectively. U >> V in this regime. In 
order map (|l]) to (||), we demand that U » J, where 
J is the Josephson coupling between islands. The easy- 
plane regime is given by V < 4J. As emphasized earlier, 
this inequality is crucial for eliminating possible charge 
density states, which may also exhibit reduced critical 
currents at commensurate fillings. 

The junction parameters reported in Ref. Q were 
Co =0.64fF , and C =1.0fF, and J =0.63°K. To reach the 
desired regime, these parameters should be modified to 
approximately C =3.0fF, and C =0.3£F, and J =0.2°K. 

Acknowlegements: Useful discussions with J. Avron, 
A. Stern, A. van Oudenaarden, and S-C. Zhang are grate- 
fully acknowledged. The authors thank the Department 
of Physics, Stanford University, where part of this work 
was performed. AA acknowledges a grant from the Israel 
Science Foundation. 



[10 
[11 

[12 

[is; 

[14 

[is; 

[16 
[17. 
[18. 

[19 

[20 

[21 
[22 



569 (1956); K.S. Liu and M.E. Fisher, J. Low. Temp. 
Phys. 10, 655 (1973); A. Aharony and A. Auerbach Phys. 
Rev. Lett. 70, 1874 (1993). 

M. Oshikawa, M. Yamanaka, I. Affleck, Phys. Rev. Lett. 
78, 1984 (1997). 

While (|l|) is sometimes used as a coarse grained theory 
of homogenous superconductors, it does not contain the 
fermion (pair breaking) excitations. 

For bipartite lattices, we replace the negative Bose hop- 
ping by a positive xy exchange using a sublattice rotation 
@. 

A. Auerbach, "Interacting Electrons and Quantum Mag- 
netism", (S pringer- Verlag. NY ) , Ch. 3; Chia Laguna lec- 
ture notes, cond-mat/9801294 . 

R. Tao and F.D.M. Haldane, Phys. Rev. B33, 3844 
(1986). 

Y. Gefen and D. J. Thouless, Phys. Rev. B47, 10423 
(1993). 

A. Auerbach, Phys. Rev. Lett. 80, 817 (1998). 
F.D.M. Haldane, Phys. Lett. A93, 464 (1983). 

D. V. Khveshchenko, Phys. Rev. B 50, 380 (1994); E. 
Dagotto and T.M. Rice, Science 271, 618 (1996); G. 
Sierra, Jour. Phys. A 29, 3299 (1996). 

The minigaps of Fig. [j] are between Sl ot = states, while 
magnons of S'i ot = ±1, have lower excitation energies due 
to easy plane anisotropy J z < J. 

E. Lieb, T. Schultz, D. Mattis, Ann. Phys. (NY.) 16, 407 
(1961). I. Affleck and E.H. Lieb, Lett. Math. Phys. 12, 
57 (1986); I. Affleck, Phys. Rev. B, 37, 5186 (1988). 
Higher order corrections were obtained for the Heisenberg 
case by S. Chakravarty, Phys. Rev. Lett., 77 4446 (1996). 
L.L. Sohn and M. Octavio, Phys. Rev. B 49, 9236 (1994); 



[1] S. Doniach, Phys. Rev. B24, 5063 (1981). 

[2] D.B. Haviland, Y. Liu and A.M. Goldman, Phys. Rev. 
Lett. 62, 2180 (1989); A.F. Hebard and M.A. Paalanen, 
Phys. Rev. Lett. 65, 927 (1990); A. Yazdani and A. Ka- 
pitulnik, Phys. Rev. Lett. 34, 3037 (1995). 

[3] P. Delsing, CD. Chen, D.B. Haviland, Y. Harada and 
T.Claeson, Phys. Rev. B50, 3959 (1994). 

[4] A. van Oudenaarden and J.E. Mooij, Phys. Rev. Lett. 76, 
4947 (1996); A. van Oudenaarden, Ph.D. Thesis, (Delft 
University Press, 1998) 

[5] M.P.A. Fisher, P.B. Weichman, G. Grinstein and D.S. 
Fisher, Phys. Rev. B40, 546 (89); C. Bruder, R. Fazio, 
A. Kampf, A. van Otterlo and G. Shon, Phys. Scr. 42, 
159 (1992); R.T. Scalettar, G.G. Batrouni, AP. Kampf 
and G.T. Zimanyi, Phys. Rev. B51, 8467 (95). 

[6] G. E. Volovik, JETP Lett. 15, 81 (1972); P. Ao and D. J. 
Thouless, Phys. Rev. Lett. 72, 132 (1994); D. P. Arovas 
and J. A. Freire, Phys. Rev. B 55, 1068 (1997). 

[7] M.P.A. Fisher and D.H. Lee, Phys. Rev. B39, 2756 
(1989). 

[8] A. Stern, Phys. Rev. B50, 10092 (1994). 

[9] T. Matsubara and H. Matsuda, Prog. Theor. Phys. 16 



4 



